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1  Abstract 


There  are  two  main  components  of  the  work  perfonned  during  this  project.  First,  a  finite-rate 
vibrational  energy  relaxation  model  was  developed  to  study  how  the  possible  presence  of  water 
vapor  in  the  Tunnel  9  facility  could  affect  the  flow  properties.  We  find  that  water  vapor  levels  at 
about  500  parts  per  million  (ppm)  could  cause  appreciable  relaxation  of  the  nitrogen  test  gas, 
however,  even  at  these  large  levels,  the  effect  on  the  flow  quality  is  minimal.  Secondly,  a  novel 
grid  generation  scheme  for  hypersonic  nozzle  flows  is  presented.  The  grid  is  based  on  the 
characteristic  lines  of  the  supersonic  regions  of  the  flow.  This  allows  for  grid  alignment  and 
clustering  along  shock  waves  present  in  the  flow.  Results  from  simulations  run  on  the  AEDC 
Hypervelocity  Tunnel  9  Mach  14  nozzle  give  good  agreement  with  simulations  performed  using 
a  conventional  grid.  Some  results  from  an  investigation  into  the  effect  of  wall  contour 
discontinuities  present  in  the  Mach  14  Nozzle  are  presented.  Avenues  for  future  research  into 
flow  irregularities  in  the  Mach  14  nozzle  are  outlined. 


2  Introduction 

Accurate  predictions  from  computational  fluid  dynamics  (CFD)  simulations  are  essential  to  the 
design  of  hypersonic  vehicles.  In  order  to  have  confidence  in  the  accuracy  of  CFD  simulations,  it 
is  necessary  to  have  experimental  data  to  validate  the  CFD  results,  for  example  measurements 
from  a  wind  tunnel.  Thus,  it  is  critical  to  ensure  that  the  wind  tunnel  design  is  of  the  highest 
quality,  and  any  non-uniformities  in  the  nozzle  flow  are  minimized.  It  is  technically  challenging 
to  design  a  wind  tunnel  that  can  produce  the  hypervelocity  flows  that  are  required  for  these 
purposes. 

The  Arnold  Engineering  Development  Center  (AEDC)  Hypervelocity  Tunnel  9  is  one  of  the 
nation's  premiere  hypersonic  test  facilities.  The  tunnel  is  large,  has  run  times  on  the  order  of 
several  seconds,  and  can  be  run  from  Mach  8  to  Mach  14  across  a  wide  range  of  unit  Reynolds 
numbers.  These  factors  make  Tunnel  9  unique  and  useful  for  obtaining  force  and  moment  test 
data  of  hypersonic  vehicle  designs.  However,  the  Mach  14  nozzle  produces  non-unifonnities  in 
the  test  section  flow  that  are  not  desirable  [1,2].  Calibration  runs  with  Pitot  pressure  rakes 
suggest  that  a  series  of  weak  shock  waves  cross  the  test  section,  which  appear  as  variations  in  the 
core  flow  Pitot  profile. 

A  potential  source  of  the  non-ideal  flow  could  be  the  presence  of  vibrational  nonequilibrium  in 
the  flow  of  nitrogen  in  the  AEDC  Tunnel  9  nozzle.  Standard  vibrational  relaxation  models 
predict  that  the  nitrogen  test  gas  vibrational  modes  freeze  close  to  the  throat.  However  careful 
CARS  (coherent  anti-Raman  scattering)  measurements  at  the  higher  operating  pressures  show 
that  there  is  no  detectable  vibrational  excitation  at  Mach  10  and  14  [3,4],  The  first  part  of  the 
research  project  studied  the  effect  of  water  vapor  on  the  vibrational  relaxation  of  the  test  gas  in 
an  attempt  to  understand  this  discrepancy  between  simulations  and  experiment. 
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Previous  numerical  studies  of  ideal  smooth-surface  flow  in  the  Mach  14  nozzle  do  not  show  the 
presence  of  weak  shock  waves  [1].  However,  it  is  known  that  the  nozzle  as-built  contour  is  not 
ideal,  and  in  fact  there  are  imperfect  joints  in  the  nozzle  that  cause  small  discontinuities  in  the 
wall  contour.  Studies  that  have  included  a  representative  model  for  similar  joint  imperfections  in 
the  Mach  8  nozzle  indicate  that  the  flow  develops  with  weak  shock  waves  passing  through  the 
test  section  [2].  Accurate  numerical  simulations  that  capture  the  effects  of  these  geometric 
imperfections  are  challenging,  as  the  waves  are  weak  as  compared  to  the  expansion  of  the  mean 
flow.  Additionally,  the  waves  are  oblique  to  a  conventional  grid,  and  thus  cross  the  grid  cells  in  a 
way  that  is  very  numerically  dissipative.  In  order  to  understand  the  influence  of  the  weak  shock 
waves,  the  grid  must  be  extremely  fine.  Due  to  the  large  range  of  length  scales  present  in  the 
problem,  even  an  axisymmetric  simulation  using  a  sufficiently  line  grid  of  the  Mach  14  nozzle 
would  be  beyond  current  computational  capability. 

This  grant  helped  develop  a  novel  grid  generation  approach  to  overcome  the  shortcomings  of 
conventional  CFD  grids  with  respect  to  this  problem,  namely  that  conventional  grids  cannot 
predict  the  Pitot  pressure  distribution  in  the  test  section  of  the  Mach  14  nozzle.  Because  the 
waves  generated  by  the  joint  imperfections  propagate  along  the  characteristic  lines  of  the  flow,  a 
grid  in  which  the  cells  are  aligned  with  the  characteristics  offers  advantages  over  a  conventional 
grid.  With  this  grid,  the  weak  shock  waves  do  not  cross  the  grid  cells  in  an  oblique  manner, 
reducing  numerical  dissipation  as  compared  to  a  conventional  grid.  Additionally,  the  grid  cells 
can  be  clustered  in  the  region  of  the  shock  waves,  which  greatly  reduces  the  number  of  elements 
required  in  the  simulation,  while  maintaining  the  necessary  accuracy  to  capture  the  effect  of  the 
waves.  This  allows  the  simulations  to  be  done  at  a  reasonable  computational  cost. 

This  final  report  summarizes  progress  toward  understanding  the  irregularities  in  the  Mach  14 
nozzle.  Results  of  simulations  aimed  at  capturing  the  effect  of  known  joint  imperfections  on  the 
nozzle  flow  will  be  shown,  and  avenues  for  future  possible  research  are  presented. 


3  Statement  of  Objectives 

1 .  Study  the  effect  of  water  vapor  on  the  vibrational  energy  relaxation  of  the  test  gas. 

2.  Development  of  grid  generation  code  which  aligns  grid  elements  with  characteristic  lines 
in  the  flow  of  the  Mach  14  nozzle 

a)  Create  lattice  of  grid  points  from  intersection  of  characteristic  lines  of  ideal  geometry 
flow 

b)  Create  axisymmetric  finite  volume  mesh  from  these  grid  points 

c)  Add  representative  models  for  nozzle  joint  imperfections  in  the  finite  volume  mesh; 
refine  where  necessary  to  capture  effect  of  joints  on  flow  physics 

3.  Simulate  flow  in  Mach  14  nozzle  with  imperfect  geometry  at  Re  =  0.5M/ft  using  US3D 

a)  Capture  individual  effects  of  both  the  step  joint  and  the  V-notch 

b)  Investigate  the  interaction  and  combined  effect  of  both  joints 
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4.  Analyze  the  data  from  the  simulations 

a)  Compare  pitot  pressure  measurements  to  simulation  results  at  various  locations  in  the 
test  cell 


4  Research  Effort 

4.1  Introduction 

The  AEDC  Hypervelocity  Wind  Tunnel  Number  9  is  located  at  the  White  Oak,  Maryland  test 
facility.  It  is  a  blowdown  tunnel  which  uses  nitrogen  as  the  working  fluid.  Tunnel  9  is  capable  of 
being  run  with  nozzles  that  can  generate  Mach  7,  8,  10,  and  14  flows  over  a  wide  range  of  unit 
Reynolds  numbers.  The  Mach  14  nozzle  has  been  the  nozzle  of  interest  for  this  research  due  to 
the  non-uniformity  present  in  the  test  cell  flow.  The  Mach  14  nozzle  is  40  feet  long  with  a  60-in. 
diameter  exit,  feeding  into  a  test  section  that  is  over  12  feet  long.  The  high  Mach  and  Reynolds 
numbers  capable  in  the  Mach  14  nozzle,  combined  with  its  large  size  and  long  run  times,  make  it 
possible  to  collect  force  and  moment  data  over  a  wide  range  of  angles  during  a  single  run.  Thus, 
making  progress  toward  understanding  the  non-uniformities  present  in  the  flow  has  been  of 
utmost  importance  for  this  research. 

Research  was  completed  in  simulating  the  Mach  14  nozzle  with  the  known  joint  imperfections 
using  the  characteristics-based  grid,  in  an  attempt  to  resolve  the  weak  shock  waves  present  in  the 
flow  and  detennine  quantitatively  their  effect  on  the  test  section.  The  simulations  were 
axisymmetric  and  performed  using  the  University  of  Minnesota's  US3D  code  [5-7],  a  finite- 
volume  Navier-Stokes  solver.  As  the  grid  generation  method  is  the  novel  aspect  of  this  research, 
the  specifics  of  the  development  of  US3D  are  not  presented  here. 

4.2  Effect  of  Water  Vapor  on  Vibrational  Relaxation 

Meador,  Townsend,  and  Miner  [8]  developed  a  model  for  the  de-excitation  of  nitrogen  in  the 
presence  of  water  vapor.  They  point  out  that  water  molecules  have  a  special  property  that  makes 
them  very  efficient  quenchers  of  vibrational  excitation.  This  is  because  of  the  geometry  of  the 
water  molecule:  two  light  hydrogen  atoms  bonded  to  a  massive  oxygen  atom  at  nearly  right 
angles.  This  makes  the  molecule’s  moment  of  inertia  small,  with  many  closely  spaced  rotational 
levels.  Thus,  there  are  many  near  resonant  energy  matches  between  the  vibrational  levels  of  most 
molecules  and  the  vibrational-rotational  levels  of  water.  In  a  hypersonic  flow  of  nitrogen  with  a 
trace  level  of  water  vapor,  the  vibrationally  excited  nitrogen  transfers  its  vibrational  energy  to  the 
vibrational-rotational  energy  of  water  vapor  through  a  very  rapid  process  (a  V-VR  transition). 
Then  the  excited  water  vapor  de-excites  through  vibrational-translational  (V-T)  transitions  due  to 
collisions  with  N2.  The  two  critical  reactions  are: 
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TV*  +H20-*N2  +H20 * 

h2o*  +n2  ->  h2o  +  n2 

Here,  the  *  superscript  indicates  a  molecule  that  is  excited  to  the  first  vibrational  state  above  the 
ground  state;  the  species  without  an  asterisk  are  ground  state  molecules.  Conventional  CFD 
models  of  vibrational  relaxation  usually  only  include  vibrational-translational  (V-T)  energy 
transitions.  They  typically  do  not  include  the  very  important  V-VR  process  shown  above,  and  as 
a  result,  under-predict  the  importance  of  water  vapor  on  vibrational  quenching. 

We  use  an  existing  computational  fluid  dynamics  code  to  study  the  effects  of  water  vapor  on  the 
Tunnel  9  Mach  14  flow  [9].  This  code  was  designed  for  the  analysis  of  hypersonic  converging- 
diverging  nozzles,  and  uses  an  excluded  volume  equation  of  state  to  model  high-pressure  non¬ 
ideal  gas  effects,  along  with  a  finite-rate  thermo-chemical  model  for  high-temperature  reacting 
flows.  As  discussed  above,  the  concentration  of  water  vapor  is  known  to  strongly  influence  the 
vibrational  relaxation  of  molecules,  including  that  of  nitrogen,  which  is  used  as  a  test  gas  in  the 
facility.  Therefore,  the  nozzle  simulation  code  was  modified  to  incorporate  the  presence  of  water 
and  its  interaction  with  the  nitrogen  test  gas.  We  use  the  model  due  to  Meador  et  al8  with  four 
chemical  species  (N2,  H20,  N2  ,  and  H20  )  and  the  following  five  reactions: 


TV*  +N2  **  N2  +  N2  (Rl) 

TV*  +  H20  TV2  +  H20  (R2) 

TV*  +H20  N2+H20*  (R3) 

H20 *  +  H20  **  H20  +  H20  (R4) 

H20*  +  N2  **  H20  +  N2  (R5) 


The  rates  of  these  reactions  were  obtained  from  three  separate  sources.  The  rate  constant  for  the 
first  reaction  (Rl)  originates  from  Vincenti  and  Kruger  [10],  and  it  can  be  derived  from  the 
relationship  between  the  relaxation  time  and  the  rate  constant  for  V-T  transitions.  The  relaxation 
time,  x,  is  the  time  required  for  the  difference  in  the  total  vibration  energy  to  the  equilibrium 
energy  to  reach  Me  of  its  initial  value.  For  a  simple  harmonic  oscillator,  we  have  the  following 
relationship  between  the  relaxation  time  and  the  rate  of  V-T  transitions  from  the  v  =  1  state  to  the 
v  =  0  state,  k\$: 


x  = 


-ejr 
-e  v 


)f 


(1) 
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Here  6V  is  the  characteristic  temperature  of  vibration  for  the  employed  species.  The  data  for 
nitrogen  and  water  are  known,  and  are  given  as:  6V  =  3395  K  for  N2  and  dv  =  2295  K  for  H2O. 

The  theoretical  relaxation  time  as  a  function  of  temperature  and  pressure  was  first  derived  by 
Landau  and  Teller,  which  at  sufficiently  low  temperatures  can  be  represented  to  good 
approximation  by: 


pr  =  Cexp((if2/r)1/3]  (2) 

where  C  and  K2  are  constants,  the  latter  of  which  depends  on  the  physical  properties  of  the 
molecule  in  question.  Constants  for  various  species  and  heat-bath  molecules  are  known,  and  for 
N2  they  are:  C  =  7.12xl0'3  atm-psec  and  K2  =  1.91xl06K. 

Introducing  the  ideal  gas  law  into  Eqs.  (1)  and  (2)  gives  the  rate  for  the  first  reaction  in  the  form 
usually  used  to  quantify  the  rate  of  vibrational  excitation  or  de-excitation.  This  expression  has 
units  of  inverse  seconds  per  atmosphere  of  pressure  of  the  relevant  species.  As  we  will  see,  we 
must  convert  these  units  to  those  required  for  the  chemical  kinetics  model.  This  reaction  rate  is 
(in  s’1  atm'1): 


km  - 


exp  (-(AT2/r)'y) 


Bass  [11]  gives  the  rate  constants  for  the  three  reactions  R2,  R4,  and  R5  as: 


kR2  =  9.0  x  10  e 


6  -28.4/7”' 


kR4  =  38.5  em'5IT'n  for  T<  410  K 


kR4  =  1 .4  x  108  for  r>  410  K 
kR5  =  2.8x  109  e-5hVT'n 

Finally,  Nagel  and  Rogovin  [12]  give  the  rate  constant  for  the  remaining  reaction  in  the  model 
(R3): 


(3) 

(4) 

(5) 

(6) 

(7) 


kR3  =  2.557  x  107  e  27IT'n  for  T<  750  K  (8) 

kR3  =  2.878  xlO9  e~69-9,r'n  f0r  t  >  750  K  (9) 
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Now  we  must  convert  these  rates  to  the  conventional  form  used  for  chemical  kinetics.  To  do  so, 
consider  a  reaction  of  the  form: 


TV*  +  S  -»•  TV2  +  S  (10) 

where  .S  is  a  collision  partner.  In  the  forward  direction  the  rate  of  formation  of  species  TV?  is 
conventionally  written  in  the  form  (with  the  reaction  rate  in  units  of  m3  /  kmole  sec). 


dp  n2 
dt 


^ 'Nl  Ps  1 
M  ,  Ms  f 

N2  ^ 


(11) 


However,  the  reaction  rates  given  above  assume  a  form  of  the  reaction  given  by: 


dNN, 

dt 


=  NJ< 


(12) 


Where  for  this  reaction,  we  have: 


k=pN,  A  exp (B  /  TV})  =  p  ,  ~^—T  A  exp (B  /  TV3) 

2  2  M  * 

n2 


(13) 


Comparing  these  expressions  shows  that  the  kinetic  rate  constant  is  given  by: 

kf  =  RT  k  (14) 

Where  the  universal  gas  constant  has  the  appropriate  units  ( R  =  0.082056  m3  atm  /  kmole  K). 

Now  we  must  detennine  the  backward  reaction  rates  for  each  of  the  reactions  given  above.  We 
can  use  the  equilibrium  constant,  which  is  given  as: 


K 


eq 


(15) 


where  kb  is  the  backward  reaction  rate.  The  equilibrium  constant  can  be  evaluated  through  the 
ratio  of  molecules  that  occupy  some  energy  state  of  vibration,  denoted  by  i,  to  the  total  number 
of  molecules  present: 


TV 


-eJkT  /  V*  -eJkT 

e  '2jF 


(16) 


where  TV)  is  the  total  population  of  energy  state  i,  TV  is  the  total  number  of  molecules,  £,  is  the 
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vibrational  energy  of  state  i,  and  k  is  the  Boltzmann  constant.  It  can  be  shown  that  the 
equilibrium  constant  for  each  reaction  is  written  as  follows: 


(Keq)m=(Keq)R2=eB^,T  (17) 

(Keq)R  3=e(^-^)/r  (18) 


(Keq)R4=(Keq)R5=ee'“>olT  (19) 

With  these  expressions,  we  can  develop  the  complete  set  of  forward  and  backward  reaction  rates 
for  each  reaction.  These  have  been  implemented  in  the  nozzle  code,  [9]  and  verified  by 
comparison  to  zero-dimensional  relaxation  results  using  the  same  kinetics  model. 

With  this  thermo-chemical  model  it  is  necessary  to  compute  the  mass  fraction  of  the  gas  in  its 
initial  state  in  the  reservoir  of  the  tunnel.  Using  the  expressions  for  the  equilibrium  constant 
leads  to  relationships  between  the  mass  fractions  of  the  ground  state  and  excited  molecules  for 
each  species: 


N2  _  H2Q  e~0,H2o/T 

cn2  ch2o 


(20) 


where  cs  is  the  mass  fraction  of  species  s.  Therefore,  given  an  assumed  mass  fraction  of  water  in 
the  Tunnel  9  test  gas  we  have  the  means  with  which  to  detennine  the  remaining  initial 
constituents  in  the  gas.  From  this  point  the  relevant  analysis  has  been  completed  to  allow  full 
implementation  of  the  chemistry  model  in  the  nozzle  simulation  code. 

The  Mach  14  Tunnel  9  nozzle  geometry  was  represented  with  a  300  by  120  computational  grid, 
including  the  portion  of  the  nozzle  upstream  of  the  throat.  The  baseline  conditions  considered 
represent  the  low  Reynolds  number  Tunnel  9  operating  conditions  of  2025  psi  (14.0  MPa)  and 
2777  R  (1543  K).  Cases  with  pure  N2  and  increasing  levels  of  H20  were  run  to  assess  the  effects 
of  water  vapor  on  the  relaxation  of  the  gas  at  these  conditions. 

Figure  1  (upper)  plots  the  nozzle  centerline  translational-rotational  and  N2  and  H20  vibrational 
temperatures  for  the  baseline  conditions;  the  simulations  were  run  with  varying  levels  of  water 
vapor  from  pure  N2  to  a  mole  fraction  of  0.005  (0.5%).  Here  the  vibrational  temperatures  are 
computed  using  the  ratio  of  the  excited  to  ground  state  mass  fractions,  as  given  in  Eq.  (20).  Note 
that  increasing  levels  of  water  vapor  do  indeed  relax  the  vibrational  modes  of  the  gas.  Also, 
because  of  the  rapid  expansion,  the  relaxation  process  occurs  in  the  throat  region,  and  the 
relaxation  process  is  frozen  for  almost  the  entire  length  of  the  nozzle.  Table  1  presents  the 
centerline  nozzle  exit  plane  conditions  for  this  case.  These  results  show  that  increasing  levels  of 
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relaxation  increase  the  translation-rotational  temperature,  as  expected.  There  is  little  effect  on  the 
resulting  centerline  Mach  number. 

It  should  be  noted  that  the  use  of  temperature  (rather  than  energy)  to  quantify  the  vibrational 
nonequilibrium  is  somewhat  misleading.  The  vibrational  energy  is  non-linear  in  vibrational 
temperature,  as  shown  in  Fig.  2.  Thus,  an  apparent  strong  level  of  vibrational  nonequilibrium 
based  on  T,  only  translates  to  a  small  level  of  energy  storage  in  the  vibrational  energy  modes. 
Once  TvN2  is  below  about  600  K,  there  is  virtually  no  energy  stored  in  the  vibrational  energy 
modes.  Furthennore,  it  is  difficult,  if  not  impossible,  to  use  CARS  to  detect  vibrational  energy  at 
low  levels  of  vibrational  energy.  The  approach  used  by  Smith  cannot  detect  vibrational  excitation 
below  about  700  K  [3,4].  Thus,  even  quite  small  levels  of  water  vapor  (0.05%  mole  fraction  or 
500  ppm)  could  be  responsible  for  the  apparent  lack  of  vibrational  nonequilibrium  in  the 
measurements  of  Smith. 

Figure  1  (lower)  plots  the  same  results  for  a  higher  stagnation  pressure  condition  of  8750  psi 
(60.3  MPa)  and  2960  R  (1644  K).  At  this  higher  pressure  condition,  the  pure  nitrogen  gas  relaxes 
more  quickly  and  lower  levels  of  water  are  required  to  drive  the  N2  vibrational  energy  close  to 
equilibrium.  Again,  small  levels  of  water  vapor  could  cause  vibrational  relaxation  of  the  Nj  gas. 

Figure  3  plots  the  computed  exit-plane  Pitot  pressure  non-dimensionalized  by  the  reservoir 
pressure  for  the  two  cases.  Note  that  the  enhanced  vibrational  relaxation  caused  by  the  presence 
of  water  vapor  does  not  affect  the  Pitot  pressure  to  any  appreciable  degree.  Thus,  if  there  is  water 
vapor  present  in  the  supply  gas,  it  does  not  affect  the  operation  of  the  facility.  This  can  be  further 
confirmed  by  considering  how  the  exit-plane  Mach  number  an  velocity  vary,  as  given  in  Tables  1 
and  2.  Therefore,  this  analysis  shows  that  the  issues  with  flow  uniformity  are  unlikely  to  be 
caused  by  vibrational  relaxation  effects. 
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Figure  1:  Centerline  translational-rotational  temperature  (left),  No  vibrational  temperature  (center) 
and  H20  vibrational  temperature  (right).  Upper:  Case  1  conditions  (2025  psi  and  2777  R);  Lower: 
Case  2  conditions  (8750  psi  and  2960  R). 


Figure  2:  N2  vibrational  energy  as 
a  function  of  temperature. 
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Figure  3.  Non-dimensionalized  Pitot  pressure  at  the 
nozzle  exit  plane  for  the  two  cases  computed;  Case  1 
(left).  Case  2  (right). 
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Table  1.  Exit-plane  conditions  for  Case  1. 


H20  Mole  Fraction 

T(  K) 

TvN2  (K) 

T,H20  (K) 

u  (m/si 

M 

0 

41.9 

1542 

— 

1784 

13.52 

0.0001 

41.1 

1213 

598 

1785 

13.52 

0.0005 

42.1 

917 

542 

1787 

13.51 

0.001 

42.3 

782 

522 

1790 

13.50 

0.005 

43.3 

534 

501 

1804 

13.43 

Table  2.  Exit-plane  conditions  for  Case  2. 


H20  Mole  Fraction 

7TK) 

T,m  (K) 

TVH20  (K) 

u  (m/s> 

M 

0 

43.3 

1047 

— 

1884 

14.04 

0.0001 

43.4 

908 

542 

1885 

14.03 

0.0005 

43.7 

679 

511 

1889 

14.01 

0.001 

44.0 

577 

503 

1893 

13.99 

0.005 

45.5 

500 

500 

1911 

13.88 

4.3  Generation  of  the  characteristics-based  grid 

The  generation  of  a  grid  aligned  with  the  characteristics  of  the  flow  is  extremely  important.  The 
unique  aspects  of  this  grid  allow  the  shocks  to  be  resolved  by  the  flow  solver  in  a  manageable 
computation  time.  In  order  to  create  the  grid,  knowledge  of  certain  flow  field  variables  is 
necessary.  To  obtain  these  variables,  the  solver  is  first  run  on  an  ideal  nozzle  geometry.  This 
provides  distributions  of  Mach  number  and  flow  angle  over  physical  space,  from  which  the 
characteristics  of  the  flow  can  be  obtained.  Note  that  by  adding  contour  discontinuities  to  the 
nozzle  wall,  the  characteristics  are  not  appreciably  affected.  The  effect  of  the  shock  waves 
generated  at  the  joints  is  much  weaker  than  the  expansion  of  the  mean  flow,  and  thus  the 
characteristics  are  dominated  by  the  nominal  nozzle  geometry. 

Once  a  solution  for  the  ideal  flow  has  been  obtained,  the  characteristic  grid  can  be  generated.  To 
begin  with,  an  initial  line  of  grid  points  perpendicular  to  the  nozzle  axis  is  created  at  a  location 
downstream  of  the  nozzle  throat.  This  is  necessary  because  the  characteristics  are  only  real  for 
regions  where  the  flow  is  locally  supersonic.  Upstream  of  this  point,  a  body-fitted  Cartesian  grid 
is  used.  Characteristic  lines  from  each  of  these  initial  grid  points  can  be  created  as  rays 
originating  from  the  point.  Each  point  has  a  left  and  right-leaning  characteristic  line.  The  angle 
of  the  characteristics  are  calculated  as  follows: 

a  =  9  ±  p  (21) 

where  a  is  the  characteristic  angle,  9  is  the  flow  angle,  and  p  is  the  Mach  angle.  The  positive  sign 
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corresponds  to  the  right  leaning  characteristic  while  the  negative  corresponds  to  the  left  leaning 
characteristic.  The  flow  angle  is  determined  by  the  ratio  of  the  local  x  and  y  components  of  fluid 
velocity,  while  the  Mach  angle  is  determined  by  the  local  Mach  number: 


d  =  arctan 


/i  =  arcsin 


(22) 

(23) 


The  slope  of  the  lines  emanating  from  a  grid  point  are  found  by  taking  the  inverse  tangent  of  the 
characteristic  angle.  The  right  leaning  line  has  a  positive  slope  and  the  left  leaning  line  a  negative 
slope.  Thus,  from  the  initial  distribution  of  grid  points  a  new  set  of  grid  points  further 
downstream  is  found  by  finding  the  intersection  of  the  calculated  characteristic  lines.  This 
process  is  repeated  marching  downstream  along  the  nozzle  axis  and  through  the  test  section.  The 
resulting  set  of  grid  points  can  then  be  used  to  create  an  unstructured  grid  of  almost  exclusively 
hexahedral  elements,  the  exception  being  elements  located  on  the  centerline  which  are 
necessarily  prismatic. 

Near  the  wall,  as  the  Mach  number  decreases,  the  characteristic  lines  curve  toward  the  wall  and 
eventually  intersect  the  wall  at  an  angle  nonnal  to  the  local  contour.  To  model  this  behavior,  the 
characteristic  lines  of  the  grid  do  not  reflect  off  of  the  physical  wall,  but  rather  somewhere 
located  in  the  Mach  number  boundary  layer.  Between  this  reflection  line  and  the  wall,  a  layer  of 
structured  hexahedral  elements  is  used  with  clustering  at  the  surface.  Because  it  is  not  possible  to 
know  exactly  where  the  reflection  line  should  be  set  ahead  of  time,  the  simulation  may  have  to 
be  run  several  times  in  order  to  ensure  that  the  characteristic  grid  elements  are  aligned  with  the 
waves  present  in  the  flow  for  the  entire  length  of  the  nozzle. 

Note  that  when  selecting  the  initial  distribution  of  grid  points,  regular  spacing  is  not  necessary. 
Clustering  of  these  initial  points  results  in  grid  cell  clustering,  and  this  can  be  done  in  such  a  way 
as  to  refine  the  characteristic  grid  at  the  location  of  the  weak  shock  waves.  As  the  characteristic 
lines  do  not  dissipate,  this  method  allows  for  grid  refinement  near  the  waves  for  the  entire  length 
of  the  nozzle  without  requiring  an  inordinate  number  of  grid  elements.  Figures  4-6  show  several 
images  of  the  characteristics-based  grid. 
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Figure  4.  Characteristic  lines  are  traced  out  once  the  flow  becomes  supersonic. 


Figure  5.  Reflection  of  characteristic  lines  at  the  nozzle  wall. 
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Figure  6.  Clustering  of  grid  elements  along  characteristic  lines. 


Simulations  have  been  conducted  to  compare  the  perfonnance  of  the  characteristics  based  grid 
with  a  conventional  grid.  The  boundary  conditions  used  are  the  same  for  all  simulations.  At  the 
inlet,  a  subsonic  boundary  condition  based  on  a  total  enthalpy  calculation  is  applied,  and  a 
supersonic  boundary  condition  is  applied  at  the  outflow.  The  no-slip  condition  is  applied  at  the 
wall  for  flow  velocity,  and  wall  temperature  is  read  in  from  a  file  and  interpolated  along  the 
nozzle  axial  direction.  All  simulations  are  axisymmetric. 

The  gas  used  in  all  simulations  is  nitrogen.  The  nozzle  reservoir  parameters  are:  pressure  of  14.7 
MPa,  temperature  of  16.39  K  and  density  of  29.23  kg/m3.  These  reservoir  parameters  correspond 
to  a  unit  Reynolds  number  of  approximately  0.5  x  106  per  foot  in  the  test  cell  as  measured 
experimentally. 

Figure  7  compares  the  simulation  results  for  ideal  flow  using  a  conventional  body-fitted 
structured  grid  and  the  characteristics  grid.  The  agreement  between  the  two  grids  for  this  case  is 
excellent,  lending  confidence  to  the  decision  to  use  the  characteristics  grid  to  simulate  the  effects 
of  the  nozzle  joint. 
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Figure  7.  Top:  Mach  number  contours  using  a  conventional  grid.  Bottom:  Characteristic  grid. 


4.4  Results 

Simulations  have  been  perfonned  in  an  attempt  to  analyze  the  effect  of  the  nozzle  joint 
imperfections.  Due  to  high  wall  temperatures  in  the  region  of  the  throat,  the  nozzle  throat  piece 
shrinks  inward  due  to  thermal  expansion.  The  stresses  associated  with  this  expansion  can  cause  a 
discontinuity  that  manifests  itself  as  a  rearwards-facing  step  that  can  grow  in  height  to  the  order 
of  1mm,  located  71.6  mm  downstream  of  the  nozzle  throat.  Simulations  were  performed  using  a 
somewhat  typical  step  height  of  0.5  mm.  Additionally,  there  is  a  .010"  wide  V-notch  located  234 
mm  downstream  of  the  throat.  Figures  8  and  9  demonstrate  how  refining  the  characteristics 
based  grid  is  done  to  resolve  the  flow  at  the  locations  of  the  joint  imperfections. 
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Figure  8.  Grid  refinement  near  the  backward  facing  step  and  along  the  wave  generated  by  the 

step. 


Figure  9.  Grid  refinement  near  the  v-notch  and  along  the  wave  generated  by  the  notch. 
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4.4.1  Presence  of  Waves  Generated  by  Joint  Imperfections 

Figure  10  shows  the  steady  state  Mach  number  distribution  for  a  simulation  run  with  the  0.5  mm 
step.  As  predicted,  the  presence  of  the  step  generates  weak  shock  waves  that  travel  down  the 
length  of  the  nozzle.  It  has  been  suspected  that  these  waves  are  responsible  for  losses  and  non¬ 
uniformities  of  Pitot  pressure  in  the  test  cell.  Figure  1 1  shows  how  the  centerline  Mach  number 
varies  with  axial  position  for  cases  run  with  the  two  joint  imperfections  individually. 

It  is  clear  that  the  waves  generated  by  the  imperfect  nozzle  geometry  cause  significant  deviations 
from  the  nominal  flow.  Figure  12  shows  centerline  Mach  and  Pitot  pressure  data  for  a  simulation 
in  which  both  joints  were  modeled  together.  This  figure  demonstrates  that  the  effect  of  both 
joints  can  be  captured  simultaneously,  and  there  appears  to  be  some  interaction  effect  between 
the  waves,  as  the  Mach  number  does  not  behave  the  same  as  in  the  isolated  cases. 


m:  0.1  1.1  2.1  3.1  4.1  5.1  6.1  7.1  8.1  9.1  10.1  11.1  12.1  13.1 


Figure  10.  Mach  number  contours  showing  the  presence  of  weak  shock  waves  generated  by  the 

step. 
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Figure  11.  Centerline  Mach  number  for  simulations  with  non-ideal  geometries. 


Centerline  Mach/Pitot  Data,  .020'  bump  +  v-notch 


X 


Figure  12.  Centerline  Mach  and  pitot  pressure  data  for  simulation  with  both  joints. 
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4.4.2  Comparison  of  Flow  Data 


Figures  13  and  14  show  the  Pitot  pressure  distributions  at  multiple  locations  in  the  test  cell  for 
the  simulation  with  both  joints  as  compared  to  experimental  data  as  well  as  simulation  data  using 
nominal  nozzle  geometry  and  a  characteristics-based  grid. 


Axial  position;  42"  from  nozzle  end 


Figure  12.  Pitot  pressure  data  42”  into  the  test  cell. 
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These  plots  demonstrate  some  key  points.  First,  the  waves  generated  by  the  joints  cause 
significant  losses  in  the  Pitot  pressure  as  compared  to  ideal  geometry.  Interestingly,  increasing 
the  size  of  the  step  height  does  not  have  a  very  significant  effect  on  the  overall  Pitot  pressure 
loss.  It  is  possible  that  there  is  a  critical  height  for  which  a  wave  will  form,  and  increasing  the 
height  further  will  change  the  location  the  wave  settles  at  rather  than  increasing  its  strength. 

Also  note  that  while  there  is  a  combinatorial  effect  of  the  waves  generated  by  the  two  joints,  it  is 
not  additive.  The  loss  in  Pitot  pressure  for  the  combined  case  is  less  than  the  sum  of  the  losses 
for  the  individual  cases.  Additionally,  there  are  still  unexplained  discrepancies  between  the 
computed  Pitot  pressure  profile  and  the  experimental  results.  There  may  be  additional  nozzle 
imperfections  causing  losses  or  variations  in  the  Pitot  pressure  that  have  not  been  taken  into 
account  during  this  study. 
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5  Summary  and  Outlook 


A  grid  generation  scheme  has  been  developed  in  which  the  grid  elements  can  be  aligned  and 
refined  near  joint  imperfections  and  the  waves  generated  by  them  in  the  Tunnel  9  Mach  14 
nozzle.  This  technique  offers  advantages  that  have  allowed  us  to  capture  the  effect  of  waves  in 
the  flow  that  have  previously  eluded  us.  This  has  led  to  a  plausible  explanation  for  some  of  the 
loss  of  Pitot  pressure  in  the  flow. 

However,  there  is  still  significant  disagreement  between  simulation  and  experimental  for  this 
data  for  this  problem.  In  order  to  fully  understand  the  test  section  core  flow,  additional  research 
into  the  tunnel  flow  needs  to  be  done.  It  may  be  necessary  to  include  more  as-built  tunnel  details 
into  the  simulations  to  accurately  capture  the  real  flow.  Some  of  these  details  may  require  a 
multidimensional  (rather  than  axisymmetric)  approach.  Overall,  current  results  are  not  definitive, 
but  they  are  promising  regarding  the  possibility  of  understanding  the  Mach  14  nozzle  flow  in  the 
future. 
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